RECENT DEVELOPMENTS IN HEAVY FLAVOUR 

PRODUCTION 

G. KRAMER 

- - - Universitat Hamburg - II. Institut fiir Theoretische Physik 

, Luruper Chaussee 149 - 22761 Hamburg - Germany 

o 

fSJ ■ We review one-particle inclusive production of heavy-flavoured hadrons in a framework 

which resums the large coUinear logarithms through the evolution of the FFs and PDFs 
and retains the full dependence on the heavy-quark mass without additional theoretical 
assumptions. We focus on presenting results for the inclusive cross section for the 
' production of charmed mesons in pp collisions and the comparison with CDF data 

from the Tevatron as well as on inclusive B-meson production and comparison with 
recent CDF data. The third topic is the production of D* mesons in photoproduction 
and comparison with recent HI data from HERA. 
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1 Introduction 



■ Onc-particlc inclusive production processes provide extensive tests of perturbative quantum 

chromodynamics (pQCD). In contrast to fully inclusive processes, it is possible to study 
distributions in the momentum of the final-state particle and to apply kincmatical cuts 
^ ' close to the experimental situation. On the other hand, contrary to even more exclusive 

'nI" , cases, QCD factorisation theorems [1, 2] still hold stating that this class of observables 

can be computed as convolutions of universal parton distribution functions (PDFs) and 
fragmentation functions (FFs) with perturbatively calculable hard scattering cross sections. 
As is well-knovifn, it is due to the factorisation property that the parton model of QCD 
has predictive power. Hence, tests of the universality of the PDFs and FFs are of crucial 
' importance for validating this QCD framework. At the same time, lowest-order expressions 

. for the hard scattering cross sections are often not sufficient for meaningful tests and the 

use of higher order computations is needed. 

The perturbative analysis is becoming more involved and interesting if the observed final 
state hadron contains a heavy (charm or bottom) quark. In this case, the heavy-quark mass 
m enters as an additional scale. Clearly, the conventional massless formalism, also known 
as zero-mass variable-flavour- number scheme (ZM-VFNS), can also be applied to this case, 
provided the hard scale Q of the process is much bigger than the heavy-quark mass so that 
terms m/Q are negligible. However, at present collider energies, most of the experimental 
data lie in the kinematic region Q ^ m and it is necessary to take the power-like mass terms 
into account in a consistent framework. 

The conventional calculational scheme is the so-called massive scheme or fixed-flavour- 
number scheme (FFNS) [3], in which the number of active flavours in the initial state is 
limited to = 3 (n^ = 4) in the case of massive charm (bottom) production, and the c 
(6) quark appears only in the final state. In this case, the c (&) quark is always treated as 
a heavy particle, not as a parton. The actual mass parameter m is explicitly taken into 
account along with px- In this scheme, m acts as a cutoff for the initial- and final sate 
coUinear singularities and sets the scale for the perturbative calculations. A factorisation 
of these would-be initial- and final state coUinear singularities is not necessary, neither is 
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the introduction of a FF for the transition b ^ B. However at NLO, terms proportional to 
Ug hi(ji^/m^), where as is the strong couphng constant, arise from coUinear gluon emissions 
by c (b) quarks or from branchings of gluons into coUinear cc (bb) pairs. These terms are 
of order 0(1) for large pr and spoil the convergence of the perturbation series. The FFNS 
with nj = 3 (n/ =4) should be limited to a rather limited range of pr, from pT = to 
Pt ^ m. The advantage of this scheme, however, is that the /p^ power terms are fully 
taken into account. 

The ZM-VFNS and FFNS are valid in complementary regions of pT, and it is desirable to 
combine them in a unified approach that incorporates the advantages of both schemes, i.e. 
to resum the large logarithms, retain the full finite-m effects, and preserve the universality 
of the FFs. An earlier approach to implement such an interpolation is the so-called fixed- 
order-next-to-leading logarithm (FONLL) scheme, in which the conventional cross section 
in the FFNS is linearly combined wit a suitably modified cross section in the ZM-VFNS 
with perturbative FFs, using a pr-dependent weight function [4]. Then the FONLL cross 
section is convoluted with a non-pertubative FF for the b ^ B transition. 

The subject of this review is the theoretical description of one-particle inclusive pro- 
duction of heavy-flavoured hadrons Xh = D, B^K,., ... in a massive variablc-flavour-number 
scheme (GM-VFNS). In such a scheme the large coUinear logarithms of the heavy-quark 
mass In ^/m are subtracted from the hard scattering cross sections and resummcd through 
the evolution of the FFs and PDFs. At the same time, finite non-logarithmic mass terms 
m/Q arc retained in the hard part and fully taken into account. 

In order to test the pQCD formalism, in particular the universality of the FFs, it is 
important to provide a description of all relevant processes in a coherent framework. There- 
fore, it is important to work out the GM-VFNS at next-to-leading order (NLO) of QCD 
for all the relevant processes. Previously, the GM-VFNS has been applied to the following 
processes: 74-7^ D*+ -I- X (direct part) [5], 7 + 7 ^ -I- X (single resolved part) [6], 
7-hp ^ D*+ +X (direct part) [7], p + p (D" , D*+ , D+ , Df) + X [8, 9, 10], where the 
latter results for hadron-hadron collisions also constitute the resolved contribution to the 
photoproduction process 7 -I- p — > X^ + X . 

In this contribution, I will review the progress achieved in describing the production 
of heavy-flavoured hadrons Xh in hadron-hadron and photon-proton collisions in the GM- 
VFN scheme as it has been worked out recently. The main focus will be on the comparison 
with experimental data from CDF at the Tevatron for p + p ^ (£)",£)*+, D+,Df) + X and 
p+p^ B+ +X and from HI at HERA for 7 + p ^ D*+ + X. 

2 Theoretical Framework 
2.1 GM-VFNS 

The differential cross sections for inclusive heavy-flavoured hadron production can be com- 
puted in the GM-VFNS according to the familiar factorisation formulae, however, with 
heavy-quark mass terms included in the hard scattering cross sections [11]. Generically, the 
physical cross sections are expressed as convolutions of PDFs for the incoming hadron(s), 
hard scattering cross sections, and FFs for the fragmentation of the outgoing partons into 
the observed hadron. All possible partonic subprocesses are taken into account. The mas- 
sive hard scattering cross sections are constructed in a way that in the limit m ^ the 
conventional ZM-VFNS is recovered. A more detailed discussion of the GM-VFNS and the 
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Figure 1: The cross section for inclusive D* production in e^e annihilation evaluated in 
NLO is compared with from CLEO [22] and BELLE [23] (left) as well as from the ALEPH 
[17] and OPAL [18] data (right). The three curves in the right figure correspond to the 
Z — > cc, Z bb and full samples. 

construction of the massive hard scattering cross sections can be found in Refs. [8, 9] and 
the conference proceedings [12, 13, 14, 15]. 

2.2 Fragmentation Functions 

A crucial ingredient entering these calculation are the non-perturbative FFs for the transition 
of the final state parton into the observed hadron X^. For charm-fiavoured mesons, Xc, such 
sets of FFs have been constructed quite some time ago. For Xc = D*^, FFs were extracted 
at LO and NLO in the MS factorisation scheme with jif ~ 5 massless quark fiavours [16] from 
the scaled-energy (x) distribution da/dx of the cross section of e+e^ D*^ + X measured 
by the ALEPH [17] and OPAL [18] collaborations at CERN LEPl. Recently, this analysis 
was extended [19] to include Xc = D'^ , , , A'^ by exploiting appropriate OPAL data 
[20]. In Refs. [16, 19], the starting scales /xg for the DGLAP evolution of the a ^ Xc FFs 
in the factorisation scale ii'p have been taken to be /io = 2mc for a = g,u,u,d,d,s,s,c,c 
and /io = 2mf, for a = b,b. The FFs for a = g,u,u,d,d,s,s were assumed to be zero 
at fip = fio and were generated through the DGLAP evolution to larger values of ^'p. 
For consistency with the MS prescription for PDFs, we repeated the fits of the Xc FFs 
for the choice = mc,mb [21]. This changes the c-quark FFs only marginally, but has 
an appreciable effect on the giuon FF, which is important at Tevatron energies, as was 
found for Z?*^ production in Ref. [8]. In the meantime much more accurate data for the 
inclusive production of D'^^D^ and D*^ mesons in e'^e" annihilation have been published 
by the CLEO [22] and the BELLE [23] collaborations. With these data new FFs have been 
constructed. These fits were done in the framework of the GM-VFNS, where the finite 
charm- and bottom-quark masses were kept in the hard scattering cross sections. A global 
fit with the data from CLEO and BELLE at 10.52 GeV together with the ALEPH [17] and 
OPAL [18] data at the Z-resonance are shown in Fig. 1 [24]. 

Already many years ago we made an analysis towards FFs for bottom-flavoured mesons 
[25] by using data from the OPAL collaboration at LEPl [26]. In the last years much 
more accurate measurements of the inclusive B meson production at the Z-resonance have 
been done by the ALEPH [27], SLD [28] and the OPAL [29] collaborations. With the data 
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Figure 2: Comparison of the ALEPH [27] (circles), SLD [28] (triangles) and OPAL [29] 
(squares) data with the NLO fits using the power ansatz. The initial factorisation scale for 
all partons is /io=4.5 GeV. 

in these references combined we have performed a new fit to obtain the FFs for q,g,b ^ B, 
where, in this case, q arc the light quarks including c. To be consistent with the starting 
scale of the PDFs the FFs of these light partons are assumed to vanish at the starting scale 
Ho = rrib = 4.5 GeV and only the b —f B FF is parametrised by the usual power ansatz at 
the starting scale fio- The FFs of the light quarks and the gluon are generated via DGLAP 
evolution at higher scales. The result of the combined fit is seen in Fig. 2. All three data 
sets are consistent with each other and the fit describes the data quite well in the whole x 
range, except possibly at rather small x [30]. 

2.3 Input Parameters 

For the numerical results presented below we have chosen the following input. For the proton 
PDFs we have employed the CTEQ6.1M PDFs from the CTEQ collaboration [31, 32] and for 
the charmed meson fragmentation functions the sets from [21]. We have set rric = 1.5 GeV, 
TOb = 5 GeV (in the case of charmed meson production), mi, — 4.5 GcV (in the case of B- 
meson production) and have used the two-loop formula for ai"^\Hfi) in the MS scheme with 
Gis^\'rnz) = 0.118. The theoretical predictions depend on three scales, the renormalisation 
scale fiR, and the initial- and final-state factorisation scales fip and fi'p, respectively. Our 
default choice for hadro- and photoproduction has been = fip = ji'p = tjit, where 
rriT = ^Jp^^+rn? is the transverse mass. Scale changes are controlled by and ^f, where 
Ci?. = fJ-R/rriT, S,F = tip/mT and ^'p = iip/mT- 

3 Hadroproduction 

A few years ago the CDF collaboration has published first cross section data for the in- 
clusive production of and mesons in pp collisions [33] obtained in Run 
II at the Tevatron at center-of-mass energies of = 1.96 TeV. The data come as distri- 
butions da/dpT with y integrated over the range \y\ < 1 and the particle and antiparticle 
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Figure 3: Comparison of the CDF data [33] with our NLO predictions for The sohd 

hne represents our default prediction obtained with fif^ = = ^'p = mx, while the dashed 
lines indicate the scale uncertainty estimated by varying fiR, ^p, and /i^ independently 
within a factor of 2 up and down relative to the central values. The right figure shows the 
data-over-thcory representation with respect to our default prediction. 



contributions are averaged. 

Our theoretical predictions in the GM-VFNS are compared with the CDF data for D* 
mesons on an absolute scale in Fig. 3 (left) and in the data-over-theory representation with 
respect to our default results in Fig. 3 (right). We find good agreement in the sense that 
the theoretical and experimental errors overlap, where the experimental results are gathered 
on the upper side of the theoretical error band, corresponding to a small value of /jr and 
large values of /xf and fx'p, the dependence being dominant in the upper pt range. As 
is evident from Fig. 3 (right), the central data points tend to overshoot the central QCD 
prediction by a factor of about 1.5 at the lower end of the considered px range, where the 
errors are largest, however. This factor is rapidly approaching unity as the value of pp is 
increased. The tendency of measurements of inclusive hadroproduction in Tevatron run II 
to prefer smaller renormalisation scales is familiar from single jets, which actually favour 
IJ-R = Pt/2 [34]. It will be interesting to compare these data with predictions using the most 
recently constructed fragmentation functions based on the BELLE and CLEO data shown 
above. For more details and a comparison with the data for the _D", Z?+, and Df mesons 
we refer to Ref. [10]. 

In the GM-VFNS framework we have also calculated the cross section distribution 
da/dpT of _B-meson hadroproduction. The calculations proceed analogously to the case 
of D mesons outlined in Ref. [8]. Now the heavy quark mass m is the b quark mass m^. 
The c quark belongs to the group of light quarks q, whose mass is put to zero. 

The NLO cross section consists of three classes of contributions. 
Class (i) contains all the partonic subprocesses with & h,b ^ B transition in the final state 
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Figure 4: Transverse-momentum distribution da/dpr of pp — > B+X at cm. energy \/S — 1.96 TeV 
integrated over tiie rapidity range \y\ < 1. Tlie contributions of class (i) (solid lines) and their gg- 
initiated parts (dashed lines) evaluated at NLO in the ZM-VFNS (upper lines) and the GM-VFNS 
(lower lines) are compared. 

that have only light partons {g, q, q) in the initial state, the possible pairings being gg, gq, 
gq, and qq. 

Class (ii) contains all the partonic subproccsses with h,b B transitions in the final state 
that also have & or 6 quarks in the initial state, the possible pairings being gb, gb, qb, qh, qb, 
qb and bb. 

Class (iii) contains all the partonic subprocesses with a g,q,q B transition in the final 
state. 

In the FFNS only the contributions of class (i) are included, but the full m dependence is 
retained. On the other hand, in the ZM-VFNS, the contributions of all the three classes 
are taken into account, but they are evaluated for m = 0. In the GM-VFNS, the class- (i) 
contribution of the FFNS is matched to the MS scheme through appropriate subtractions 
of would-be coUinear singularities, and is then combined with the class- (ii) and class- (iii) 
contributions of the ZM-VFNS. Thus only the hard scattering cross sections of class (i) 
carry explicit m dependence. Specifically, the subtractions affect initial states involving 
g bb splittings and final states involving g —> bb, b gb and b gb splittings, and 
they introduce logarithmic dependences on the initial- and final-state factorisation scales in 
the hard-scattering cross sections of class (i) , which are compensated through NLO by the 
respective factorisation scale dependences by the 6-quark PDF and the b ^ B FF, respec- 
tively. It turns out that the g-quark fragmentation contribution is negligible. However, the 
gluon fragmentation reaches approximately 50% at small values of px, and somewhat less 
towards larger values of pr- 

The explicit contributions to the hard scattering cross sections of class (i) as they con- 
tribute to the final result, after all the subtractions are made, are shown in Fig. 4. The 
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Figure 5: Transverse-momentum distribution da/dpr of pp B+X at cm. energy ^/5 = 1.96 TeV 
integrated over the rapidity range |y| < 1. The central NLO prediction with = = 1 (solid 
line) of the GM-VFNS is compared with CDF data from Refs. [35] (open squares) and [36] (solid 
squares). The maximum and minimum values obtained by independently varying and £,f in the 
range 1/2 < ^r,£,f < 2 with the constraint that 1/2 < ^r/(,f < 2 are also indicated (dashed lines). 



results for m = and finite m are shown in this figure as the upper and lower solid fines, 
respectively. They constitute parts of the final ZM-VFNS and GM-VFNS results. In both 
cases, the contributions of classes (ii) and (iii) for m = still must be added to obtain the 
full predictions to be compared with experimental data. The class- (i) contributions in the 
ZM-VFNS and GM-VFNS schemes are, therefore, entitled to be negative and they indeed 
are, for px ^ 76 GeV and pr ^ 10 GeV, respectively, as may be seen from Fig. 4. Comparing 
the ZM-VFNS and GM-VFNS results, we notice that the finite-m effects are significant for 
Pt ~ 10 GeV and even cause a sign change for 10 GeV ^ ~ 76 GeV. However, as will 
become apparent below, the contributions of class (i) are overwhelmed by those of classes 
(ii) and (iii), so that the finite-m effects are washed out in the final predictions, except for 
very small values of px- It is instructive to study the relative importance of the 5(7- initiated 
contributions. They are also included in Fig. 4 for m = and finite m as the upper and lower 
dashed lines, respectively. They exhibit a similar pattern as the full class- (i) contributions 
and dominate the latter in the small-pT range. Comparing Fig. 4 with Fig. 2(c) in Ref. [8], 
we observe that the relative influence of the finite-m effects is much smaller in the c-quark 
case, as expected because the c quark is much lighter than the b quark. One can also see 
from Fig. 4 that the difference of the class (i) contributions in the GM-VFNS and ZM-VFNS 
decrease with increasing pT- 

In Fig. 5 we show the comparison of the final prediction, in which all contributions of 
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classes (i), (ii) and (iii) are combined [30] with recent Tevatron data. We compare our 
prediction to the more recent CDF data from rmi II in Refs. [35] (open squares) and [36] 
(soHd squares). In this figure the sohd line presents the central prediction for ^/j = = 1 
and the dashed lines indicate the maximum and minimum values obtained by independently 
varying and = Cf ^^le range 1/2 < S^r, $f < 2 with the constraint 1/2 < ^r/£,f < 2. 
The maximum and minimum values correspond to = 2 and ^f = 1/2, respectively. The 
variation with ^r is milder than the one with ^f- For ^f < 1, M-F reaches the starting scale 
Ho = m for the DGLAP evolution of the FFs and the 6-quark PDF at pT = ™6-\/1/Cf ~ 1- 
For smaller values of pt, there is no prediction because the FFs and the &-quark PDF are 
put to zero for fiF < ^J■o■ This explains why the pt distribution for ^f = 1/2 only starts at 
Pt = V^rrib w 7.8 GeV. The most recent data [36] nicely agree with the GM-VFNS result. 
They lie close to the central prediction, with a tendency to fall below it in the lower pT range, 
and they are comfortably contained within the theoretical error band. We conclude from 
this, that the notorious Tevatron i?- meson anomaly with data-to-theory ratios of typically 
2-3, that has been in the literature for more than a decade, is actually not present thanks 
to both experimental and theoretical progress. The previous CDF data [35] based on the 
measurement of J/ip + X final states are compatible with the latest ones for < 12 GeV, 
but are systematically below them for the larger values of pT- This inconsistency becomes 
even more apparent by noticing that Fig. 4 only contains 4 out of the 13 data points for 
Pt > 12 GeV quoted in Ref. [35] and that the omitted data points line up with the selected 
ones. This suggests that the systematical errors in Ref. [35] and perhaps also in ref. [36], 
might be underestimated and that the overall normalisation might need some adjustment. 

The measured pt distributions of Ref. [35] reaches down to almost pT = and exhibits 
a maximum at pT ~ 2.5 GeV. This small-pT behaviour is correctly reproduced in the 
FFNS without DGLAP-evolved FFs, which receive only contributions of class (i) without 
any subtractions. It is clear that our present implementation of the GM-VFNS is not 
suitable for cross sections in the small-pT region. Although the GM-VFNS is designed to 
approach the FFNS in its region of validity without introducing additional matching factors, 
to implement this numerically is not easy due to necessary cancellations between different 
terms in the calculation. The problem to achieve such cancellations is complicated by the 
extra factorisation scale; to obtain a smooth transition from the GM-VFNS to the FFNS, 
one has to carefully match terms that are taken into account at fixed order with terms 
that are resummcd to higher orders in the PDFs and FFs. In addition, it remains to be 
investigated whether a proper scale choice in the small-p^ range is required and helpful to 
ensure that the FFs and fe-quark PDF are sufficiently suppressed already at pT = 0{m). 

We extend our numerical analysis to include the NLO prediction in the FFNS with Uf = A 
massless quark flavours in the initial state, which allows us to also compare with the small-p^ 
data from Ref. [35]. We evaluate ai"'' Vfl) with n/ = 4 and A^'^'> ^ 326 MeV [31], while we 
continue using the CTEQ6.1M proton PDFs [31], in want of a rigorous FFNS set with nj = A. 
In the FFNS, there is no room for DGLAP-evolved FFs, and only b,b ^ B transitions are 
included. For simplicity, we identify b (anti)quarks with B mesons and account for non- 
perturbative effects by including the branching fraction B{b B) = 39.8% [37] as an 
overall normalisation factor, i.e. we use a. b ^ B FF of the form D{x) = B{b B)S{1 — x), 
while the g,q,q ^ B FFs are put to zero. In Fig. 6, the central FFNS (dot-dashed line), 
ZM-VFNS (dashed line), and GM-VFNS (sohd fine) predictions, for Cii = = 1, are 
compared with the CDF data from Refs. [35, 36]. As in Fig. 4, some of the data points with 
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Figure 6: Transverse-momentum distribution da/dpr of pp B+X at cm. energy ^/S = 1.96 TeV 
integrated over the rapidity range \y\ < 1. The central NLO predictions in the FFNS with Uf = 4 
and without FFs (dot-dashed line), the ZM-VFNS (dashed line), and the GM-VFNS (solid line) 
are compared with CDF data from Refs. [35] (open squares) and [36] (solid squares). 



Pt > 7 GcV from Rcf. [35] are omitted for clarity. Since the ZM-VFNS and our present 
implementation of the GM-VFNS are not applicable to the small-pr range, we show the 
respective predictions only for px > 2m = 9 GeV. The GM-VFNS prediction shown in Fig. 6 
is identical with the central one in Fig. 5. By construction, it merges with the ZM-VFNS 
prediction with increasing value of pt- In accordance with the expectation expressed in the 
discussion of Fig. 4, the difference between the GM-VFNS and ZM-VFNS results is rather 
modest also at px ^ 2m, since the m-dependent contribution, of class (i), is numerically 
small and overwhelmed by the m-independent ones, of classes (ii) and (iii). The FFNS 
prediction faithfully describes the peak structure exhibited by the next-to-latest CDF data 
[35] in the small-pT range and it also nicely agrees with the latest CDF data [36] way out to 
the largest pr values. In fact, for pT > 4m, where its perturbative stability is jeopardised by 
unresummed logarithms of the form ln(m|./m^) > 3, the FFNS prediction almost coincides 
with the GM-VFNS one, where such large logarithms are resummed. This might be a pure 
coincidence, which becomes even more apparent if we also recall that the implementation 
of the b,b ^ B transition in the FFNS is not based on a factorisation theorem and quite 
inappropriate for such large values of pr ■ 

4 Photoproduction 

Inclusive photoproduction of D* mesons, 7 + p ^ D* + X, has been studied in Ref. [7] 
where the direct part has been computed in the GM-VFNS whereas the resolved part has 
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Figure 7: Inclusive D* cross sections as a function of pt{D*) (left) and rj{D*) (right) com- 
pared to NLO QCD calculations of FMNR [42] in the FFNS and GM-VFNS for photopro- 
duction in the laboratory frame. The FMNR bottom contribution is shown separately for 
the pt{D*) distribution. 



been included in the ZM-VFNS. In this analysis the FFs of Ref. [16] and, for the resolved 
contribution, the GRV92 photon PDFs [38] have been utilized. The other parameters have 
been chosen as specified in Sec. 2.3. In Fig. 6 of Ref. [7], the central numerical predictions 
for the pt distributions of the D* meson have been compared with preliminary ZEUS data 
[39]. There exist similar data by the HI collaboration [40] which have not been used in 
this analysis. As can be seen in this figure, the agreement of the pT-distributions with 
the data is quite good down to pr — 2rac and the mass effects turn out to be small. In 
order to extend the range of applicability of the GM-VFNS into the region < 3 GeV 
more work on the matching to the 3-fixed flavour theory would be needed. Figs. 7 - 9 of 
Ref. [7], showing results for the rapidity (y), invariant mass {W) and inelasticity {z{D*)) 
distributions, have to be taken with a grain of salt since they receive large contributions from 
the transverse momentum region 1.9 < < 3 GeV which is outside the range of validity 
of the present theory. With the work in Ref. [8] , it was possible to include also the resolved 
part in the GM-VFNS. This has been done and the predictions in the complete GM-VFNS 
framework at NLO, combined with updated FFs [21], have been compared with recent HI 
photoproduction data [41]. The results of the calculation and the comparison with the data 
is shown in Fig. 7. In this figure also the predictions in the FFNS based on the FMNR 
program [42] are shown. The experimental cross section as a function of pr falls steeply 
with increasing pT as predicted by both calculations. FMNR predicts a distribution which 
decreases less steeply at large pr than the data as is seen more clearly from the plot of the 
ratio of the theoretical over the measured cross section. Also in Fig. 7 the differential cross 
section as a function of the pseudorapidity Tj{D*) is shown. This cross section decreases 
with increasing i]. Both calculations predict a similar shape and agree nicely with the data. 
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The GM-VFNS prediction shows a larger scale dependence. Otherwise the two calculations 
give rather similar results, which is remarkable, considering the very different ingredients of 
the two approaches. 

5 Summary 

We have discussed one-particle inclusive production of heavy-flavoured hadrons in hadron- 
hadron and photon-proton collisions in a massive variable-flavour-number scheme (GM- 
VFNS). The importance of a unified treatment of all these processes, based on QCD factori- 
sation theorems, has been emphasised, in order to provide meaningful tests of the universality 
of the FFs and hence of QCD. At the same time, it is necessary to incorporate heavy-quark 
mass effects in the formalism since many of the present experimental data points lie in a 
kinematical region where the hard scale of the process is not much larger than the heavy- 
quark mass. This is achieved in the GM-VFNS, which includes heavy-quark mass effects 
and still relies on QCD factorization. We have discussed numerical results for two reactions. 
In general, the description of the transverse momentum spectra is quite good down to trans- 
verse momenta px — 2to. Extending the range of applicability of our scheme to smaller pT 
would require more work on the matching to the corresponding theories in the fixed flavor 
number scheme. 
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